numpy indexing

Official documents

hr.icon

Tips

Numpy配列のAdvanced Indexingは知らないとわけわからなくなる。([], compare, condition, 条件式, fancy indexing ) - Qiita

MetPy Mondays '#90 - What is a NumPy Axis? : Unidata Developer's Blog

hr.icon

Examples

fに周波数、Pにスペクトル密度が入っているときに、Pが大きい順に対応する上位10個の周波数

code:python

code:python

import numpy as np

# lon: longitude lat:latittude

xrange=np.logical_and(lon>=120.,lon<=180.)

yrange=np.logical_and(lat>=20.,lat<=45.)

range selection 2

code:python

import numpy as np

# lon: longitude lat:latittude

xrange=((lon>136.) & (lon<140.))

yrange=((lat>36) & (lat<38))

range selection 3

code:python

import numpy as np

# lon: longitude lat:latittude

find nearest value in numpy array

code:python

import numpy as np

def find_nearest(array,value):

idx=(np.abs(array-value)).argmin()

find index where 2-dimensioal (multi-dimensional) array has maximum

code:python

ans=np.where(a == a.max())

i,j = np.unravel_index(a.argmax(), a.shape)

MetPy Mondays `#48 - Satellite Data with CLASS Part 5 : Unidata Developer's Blog

配列walkが初めて10以上になるインデックス

code:python

Trueは最大値扱いである。

NOT

否定 "-" を使う

code:python

a=np.arange(5)

Tips

...

code:python

hr.icon

Useful methods

argmax()

argmin()

searchsorted()

argsort()

argwhere()